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Abstract. We study the pseudogap Bose-Fermi Anderson model with a continuous-time 
quantum Monte Carlo (CT-QMC) method. We discuss some delicate aspects of the 
transformation from this model to the Bose-Fermi Kondo model. We show that the CT-QMC 
method can be used at sufficiently low temperatures to access the quantum critical properties 
of these models. 



Over the past decade intermetallic compounds have served as model systems to study 
instabilities of metallic magnets near zero temperature [1-3]. In particular the quantum critical 
properties found in several heavy fermion compounds seem to be beyond the Ginzburg-Landau 
paradigm of criticality [4-7]. A natural explanation for the new type of criticality invokes 
inherently-quantum critical modes in addition to the gapless order parameter fluctuations. In 
the case of heavy fermion compounds, the additional critical mode has been identified with the 
destruction of the Kondo effect [8-10]. Since the destruction of the Kondo effect is local in space 
it can be systematically studied in simplified, local models. In these quantum impurity systems 
the critical state only appears at the (spatial) boundary of suitable hosts. 

In the present article, we study the Kondo-destroying quantum critical point in the pseudogap 
Bose-Fermi Kondo model (PBFKM) with Ising anisotropy, defined as 



Where e^, uo q are the fermonic and bosonic bath dispersions, J is the Kondo coupling between 
the spin of the conduction electrons and the spin of the impurity, a is a vector of Pauli spin 
matrices and g is the coupling between the z-component of the impurity spin and the bosonic 
bath. We take a pseudogap density of states (DOS) for the conduction electrons, p c (e) oc |e| r 
for |e| < D and < r < 1/2 where we have taken the Fermi energy to be zero, and a sub-ohmic 
density of states for the bosons, pb{oj) = ^ q [5{<jj — Lo q ) — 5{<jj+LOq)\ oc sgn(w)|a;| Q! up to a cutoff A. 
A perturbative renormalization group study has been used carried out for the Bose-Fermi Kondo 
model with non-zero conduction-electron DOS (r = 0) in both the Ising and continuous spin 
symmetry cases [11]. For the pseudogapped case (r ^ 0), such a study has been possible only for 
the continuous spin symmetry case [12]. The numerical renormalization group method, on the 
other hand, has been used to study the model in the Ising limit [13]. Eq. (1) has as particular 
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limits the pseudogap Kondo and the Bose- Fermi Kondo models, which have recently been studied 
using related methods [14, 15]. Despite the seeming simplicity of Eq. (1), it is hard to study 
the quantum critical properties of the PBFKM directly. In particular, the quantum-relaxational 
regime (Tito < ksT) is difficult to address theoretically. The quantum-relaxational regime of the 
pseudogap Kondo model has been studied previously [15], by applying a CT-QMC algorithm [16- 
18] to the pseudogap Anderson model whose low-energy sector in the local moment regime can be 
mapped onto the pseudogap Kondo model via a Schrieffer- Wolff (SW) transformation. Since the 
Anderson model only involves fermions, standard perturbative expansion methods can be applied 
to it. In the CT-QMC approach, the full perturbative expansion in terms of the hybridization 
between conduction and localized fermions is then sampled stochastically using a Monte Carlo 
algorithm [16-18]. The corresponding Bose- Fermi Anderson model is 
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where is the energy level of the impurity, U is the on-site interaction , n a = d\d a , V k is the 
hybridization of the impurity with the conduction electrons, and S z = ^(n^ — n^j. The CT-QMC 
method has been extended to treat ohmic bosonic baths coupled to the charge of the impurity 
by invoking a Firsov-Lang (FL) transformation [19, 20]. Here we extend this approach to treat 
a sub-ohmic bosonic bath that couples to the spin of the impurity and explore the possibility of 
using it to access the quantum critical properties. As it turns out, the generators of the FL and 
SW transformations do not commute, raising the important question as to which is the proper 
order of applying the two transformations. 

We focus on the particle-hole symmetric case, U = — 2e^. First, we perform the FL 
transformation to eliminate the term linear in (j) exactly. We choose a generator Sfl = 
gS z J2 q ~~ an d the transformed Hamiltonian H = e SpL He~ SpL is 

H = S e k c l,a°k,a + e~d(n t + nj.) + Un^ + 5^(Vfcd£c fc)(T + V k *c{ a d a ) + ^2u q <f>\<f> q , (3) 

where d\ = 4exp(f £ g - *_,)), a = ±1 for |/|, U = U + \g 2 E g and ~e d = 

Note that h a = n a and the transformation does not destroy particle-hole symmetry. 

A modified SW transformation [21] is used to eliminate the hybridization term in Eq. (3). 
We only consider the Kondo limit, and for simplicity neglect the k dependence of V k = 
V. The generator Ssw is standard but with U,ed,d,d^ replaced by U ,€d,d,d^ , namely 
Ssw = E fc , CT ^(^f + e - d +frl e J (4<V - 4,<A)- Writing the Hamiltonian in Eq. (3) as 
H = H + H b + H h where H = e d (n t + nj.) + Unfn± + J2k,a £kc\ a c k ^, H b = J2 q u q 4>\4> q and 

H h = Efe, CT ^(4cfe,a + 4,a^) ! we have H> = e s swH e Ssw ~ H +H b +[Ssw,H b ] + ±[Ssw,Hh]- 
Projecting out unoccupied and doubly occupied states we arrive at 

H' = J2 € k C k,a C k,a + E + E i\^k,k' + \h,k'Y{^k< ,o 

k,cr q k,k',<J 

- E4 fc 4(4,^" + v^) + ^^), (4) 
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where W k w = V 2 (—^ + — ^-), J k k' = V 2 ( 1—- + I—- K ^) is the 

Kondo coupling, 5+ = 5+exp( 5 £ ? ^\ ~ $-,)), S~ = S-exp(-</ *-,)), 



Figure 1. The single particle 
Green function G(t, T) and the spin 
susceptibility x(r, T) for U = g = 
0, T = O.LD and r = 0.4 at T = 
6.6 x 10~ 4 D. Comparing the exact 
calculation and the CT-QMC result 
, „ , : we see agreement for both G and %■ 
X/P x/p 

5 = \ J2a,p <$ a a a ^dp and = \ J2a,/3 c \ cPctfi c k'fi- The third term in equation (4) represents 
a potential scattering of the conduction electrons, and the last is the Kondo term, but with 
renormalized impurity spin flip operators due to the presence of the bosonic bath. Jk t k' differs 
from the standard expression in that U and is replaced by U and e^. 

We will now discuss the opposite order of transformations, namely, H" = 
e SpL e Ssw He~ Ssw e~ SpL . Applying the SW transformation, projecting out charge fluctuations 
and then applying the FL transformation arrives at equation (4), however with Wk,k'i Jk,k' 
replaced by Wk,k', Jk,k'- We see that applying first the SW transformation, which is not ex- 
act, completely ignores the bosonic baths' influence on the charge degrees of freedom of the 
impurity. Whereas applying the FL transformation first, which is exact, correctly captures the 
bosonic baths' influence on the Anderson model which lowers the Kondo coupling. The non- 
commutativity of the two transformations quantitatively affects the effective Kondo scale at the 
quantum critical point, but it does not change the universal scaling behavior of the quantum 
critical properties because the critical value of the Kondo coupling is not universal. 

As a check on the CT-QMC approach we first compare the single particle Green function 
and the spin susceptibility for the numerical result with U = g = to the analytic result. For 
g = the bosonic bath decouples from the problem and can be ignored, and taking U = 
reduces the Hamiltonian in Eq. (2) to the resonant level model with a pseudogap. The impurity 

single particle Green function, {T T d{T)S (0)) is then G(oj) = (oj — — J2k Using 
Vk = V and taking the infinite bandwidth limit we can perform the sum over k [22] . We obtain 
G(u) = (uj — — where ReS(w) = — r(cj)tan(^)sgn(w), ImS(a;) = — T(oj) and we have 

defined the dynamic hybridization function to be T(lj) = tt\V | 2 poM r . The imaginary time Green 
function can then be obtain by Fourier transform, G(t,(3) = ^j^^— -Im(G(a; + i0 + )) and 
the local spin susceptibility can be constructed xi T i 0) = —\G{t, (3)G{—t, 0). As seen in figure 1 
we obtain quantitative agreement, within numerical accuracy, in the long time behavior for both 
G(t,0) andxW). 

We now turn to the quantum-critical properties of the PBFKM defined in Eq. (2) by 
measuring the static spin susceptibility, described below. In what follows we fix r = 0.4 and 
a = 0.6. After the (FL) transformation, the Hamiltonian in Eq. (3) can be expanded in the 
hybridization term. We use the CT-QMC algorithm to calculate the partition function, the 
single particle Green function and the local spin susceptibility. Fixing U = 0.025-D and varying 
g we can tune the model to a quantum critical point (QCP). Within the CT-QMC approach 
we measure the local spin susceptibility x(r, 0) = (T t S z (t)S z {0)} and then calculate the static 
susceptibility Xstat(/3) = Jo ^ T ^( r > 0)> where we have set the Lande g-factor and Bohr magneton 
to unity. For small g, the finite U gives rise to a Kondo screened local moment; the static 
susceptibility approaches a constant for temperatures well below the Kondo temperature, Tk- 
For large g, the impurity spin decouples from the conduction band and follows the fluctuations 
of the bosonic bath; the static susceptibility takes the Curie- Weiss form, Xstat (T) ~ T _1 . At the 
QCP, the bosonic bath acts to decohere and destroy the Kondo effect [14]. Consequently, at the 




g = 0.18D-»- 
\ f(x) = 0.42X- 0609 




Figure 2. Static susceptibility in the 
vicinity of the quantum critical point, 
g c ps 0.18D. For r = the Kondo 
temperature is, T K ps 0.029D. Well 
below Tk we see divergence of the static 
susceptibility as Xstat ~ 7"-°- 609 . 



QCP the scaling of Xstat(7") acquires an anomalous exponent Xstat(^) ~ T~ a for temperatures 
well below Tk- As seen in figure 2, using the CT-QMC approach we obtain g c ps 0.18-D and in the 
vicinity of the quantum critical point, Xstat(^) ~ T~ x with x = 0.609. Our calculated exponent 
agrees with the numerical renormalization group result within numerical accuracy [13]; the same 
exponent is also expected in related pseudogap Bose-Fermi Kondo model with continuous spin 
symmetry [12] or the Bose-Fermi Kondo model with Ising symmetry but with r = [11]. 

In conclusion, we have shown that the low energy properties of the pseudogap Bose-Fermi 
Kondo model can be addressed within a continuous-time quantum Monte Carlo approach. We 
have demonstrated that this approach correctly reproduces the exactly solvable limit of the 
pseudogap resonant level model, and been able to determine the critical behavior of the static 
local spin susceptibility in an interacting case. 
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